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Abstract 

In this paper we present an efficient computational and symbolic algorithms for solving a backward 
pentadiagonal linear systems. The implementation of the algorithms using Computer Algebra Systems 
(CAS) such as MAPLE, MACSYMA, MATHEMATICA, and MATLAB are straightforward. An ex- 
amples are given in order to illustrate the algorithms. The symbolic algorithm is competitive the other 
methods for solving a backward pentadiagonal linear systems. 
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1 Introduction 

Many problems in mathematics and applied science require the solution of linear systems having a backward 
pentadiagonal coefficient matrices. This kind of linear system arise in many fields of numerical computation 
[1, 2]. This article is a general case of the author article [3]. 



The main goal of the current paper is to develop an efficient algorithms for solving a general backward 
pentadiagonal linear systems of the form: 
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X = (xi,x 2 , ■ ■ ■ ,x n ) T , Y = (2/1,1/2, ■ ■ ■ ,Un) T and n > 5. 



A general nxn backward pentadiagonal matrix A of the form (1.2) can be stored in 5n — 6 memory locations 
by using five vectors a = (ai,a 2 , . • . ,a n _ 2 ), a = (01,03, ■ ■ . ,a„-i), b = (b 2 ,b 3 , . . .,&„), b = (63, 64, . . . ,b n ), 
and d — (di,d2, . . . ,d n ). When considering the system (1.1) it is advantageous to introduce three addi- 
tional vectors a — (a±, a.2, . . . , a n -\), f3 — {fix-, 02, ■ ■ ■ , /3 n ) an d 7 = (72, 73, • ■ • , 7n)- These vectors are 
related to the vectors a, a, d, b, and b. 



The current paper is organized as follows. In section 2, the main results are given. An illustrative 
examples are presented in section 3. In section 4, a Conclusion is given. 

'Home Address: Mansoura University, Mansoura 35516, Egypt, e-mail: abibka@mans.edu. eg 



1 



2 Main results 



In this section we are going to formulate a modified computational algorithm and an efficient symbolic 
algorithm for solving a general backward pentadiagonal linear systems of the form (1.1). To do this we 
begin by translate the system (1.1) to the following pentadiagonal linear system: 

A X X = Yi (2.1) 

where 
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Yi = (Vn, Vn-i, Vi) T and n > 5. 

Now considering the LU decomposition [4] of the matrix A\ in the form: 
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From (2.3) we obtain 
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It is not difficult to prove that the LU decomposition (2.3) exists only if fa ^ 0, i = l(l)n — 1. Moreover 
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a general backward pentadiagonal linear system (1.1) possesses a unique solution if, in addition, n ^ 0. 
On the other hand, the Determinant of the matrix Ai is given by 

n 

det(A 1 ) = '[[l3 i , (2.7) 

i=l 

and this shown the importance of the vector (3 [5] . 
We may now formulate the following results. 

Algorithm 2.1 To solve the general backward pentadiagonal linear system (1.1), we may proceed as fol- 
lows: 

step 1: Set 0i = d n , 72 = cti = K, 02 = d n -i - ai * 72, and a 2 = b n -i - 12 * K- 

step 2: If 0i= Or 2 = 0, then OUTPUT ('the method is fails'); STOP, 
step 3: For i = 3, 4, . . . , n — 1 Compute 



74 = ^ ' 

Pi ^"n—i-\-l ~~ (3 i _ 2 ^+3 ^i— lTi? 

//A = 0, tften OUTPUT('the method is fails'); STOP, 
step 4 : Compute 

0n=di- - a„_i7„. 

step 5: Set z\ = y n , z 2 = y n -i - I2Z1. 
step 6: For i = 3, 4, . . . , n Compute 

Zi = y n -i+i - a p^ 1 z -i-2 - iiZi-i. 
step 1: Compute the solution vector x using 

■r — £a. T , — z n -i-a n -ix„ 
— 0n ' n ~ 1 Pn-l 7 

For i = n — 2, n — 3, ...,1 compute 

The new algorithm 2.1 will be referred to as KBPENTA algorithm. KBPENTA algorithm for solving 
the backward pentadiagonal system (1.1) is generally preferable because the conditions 0i ^ 0, i — l(l)n 
are sufficient for the validity of it. The advantage of the vector f3 is now clear. 

The following symbolic algorithm is developed in order to remove the cases where the numeric algorithm 
KBPENTA fails. 

Algorithm 2.2 To solve the general backward pentadiagonal linear system (1.1), we may proceed as fol- 
lows: 

step 1: Set 0i = d n . 

step 2: Set 0i — x(x is just a symbolic name) whenever 0i = 0. 

step 3: Set 72 = ^j^- , «i = b n , 2 = d„_i - 0172, and a 2 = b„_i - 72&„. 

step 4 : Set 2 — x whenever 2 = 0. 

step 5: For i = 3, 4, . . . , n — 1 Compute 

On-i+l- "p7^ 1 "i-2 



7* 



ft. 



' J n-i+l — Jibn-i+2, 



0i d n —i+i ^+3 ^i — lTi 

Set 0i — x whenever t = 0. 
step 6; Compute 

7" ^ /3„_! ~ ' 

n =di- -p^h - a„_i7„. 
step 7; 5et /3„ = a; whenever n = 0. 
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step 8: Set z\ = y n , z 2 = y n -i - J2Z1. 
step 9: For i = 3, 4, . . . , n Compute 

step 10: Compute the solution vector x using 



For i = n — 2,n — 3, ...,1 compute 

1 ~ ft 

step 11: Substitute x = in all expressions of the solution vector Xi,i= 1, 2, . . . , n. 



The symbolic algorithm 2.2 will be referred to as KSBPENTA algorithm. 

In [6], Claerbout showed that the two-dimensional Laplacian operator, which appears in 3-D finite-difference 
migration, has the form of pentadiagonal matrix. If we choice di = —4, i = l(l)n, — hi — di = bi = 1, 
and yi = yj, i = l(l)n, j = n(— 1)1, we can obtain it. 



3 An Illustrative Examples 

In this section we are going to give an illustrative examples 

Example 3.1. Solve the backward pentadiagonal linear system of size 5 given by 
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(3.1) 



by using the KBPENTA algorithm and KSBPENTA algorithm. 
Solution 

(i) The application of the KBPENTA algorithm gives: 



• /3i = —1, 72 = —4, ct\ = l,p2 = 2, and a 2 = 6 (Step 1). 



. [73,74] = [2, f,-§], [/?3,/3 4 ] = [-7,f], [a 3 ,a 4 ] = [-3,f] (Step 3). 

. 75 = -§, 0b = f (Step 4). 

• z\ = 4, and z 2 = 30 (Step 5). 

• [z 3 ,z 4 ,z 5 ] = [-28,28,f ] (Step 6). 

• x 5 = 5, x 4 = 4, [xi,X2,x 3 ] = [1, 2, 3] (Step 7). 

Also the determinant of the matrix A\ is det{A\) = 160 by using (2.7). 

(ii) The application of the KSBPENTA algorithm gives: 
X := backjpenta(b, 6, d, a, a, y) = [1, 2, 3, 4, 5]. 
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Example 3.2. Solve the backward pentadiagonal linear system of size 5 given by 
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(3.2) 



by using the KBPENTA algorithm and KSBPENTA algorithm. 
Solution 

(i) The application of the KBPENTA algorithm gives: 

The method is break down since f}\ = d 5 = 0. 

(ii) The application of the KSBPENTA algorithm gives: 

X := back.penta(b, b, d, a, ~a, y)=[-^n, ^E#, ffff 



x=0 



-[1, 2, 3, 4, 5]. 

Also the determinant of the matrix A\ is det(Ai) = 88 and for more details about how to call this procedure, 
see appendix 1. 

4 Conclusion 

The methods described here are a very effective, provided optimal LU factorization is used. Our symbolic 
algorithm is competitive the other methods for solving a backward pentadiagonal linear systems which 
appear in many applications. 
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Appendix 1. A Maple procedure for solving a backward pentadiagonal linear systems 
Notes: The procedure is based on the following results: 



where 



_ Zi-aiXi+l-b n -i + ix i+2 • 
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Zi = Vn-i+i - JiZt-i - "^J* 1 Zi-2, i = 3,4:,...,n,z 1 =y n ,z 2 = y n -i - 72^1, 



H = 



-, i = 3,4,.. .,n,7 2 = 



and 



A = <i„_i + i - ^-^-&„_ i+3 - Qi_i7i, i = 3,4, .. .,n,(3i = d n ,f3 2 = d n -i - ai72, 
on = K-i+i - liK-i+2, i = 2, 3, . . . , n - 1, a\ = b n . 
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> ■#■ A Maple Procedure. 

> # Written by Dr. A. A. Karawia. 

> # To compute the solution of A general backward pentadiagonal Linear systems. 

> restart: 

> with(linalg, vector ,vectdim): 

> back_penta:=proc(bb::vector,b::vector, d::vector,a::vector,aa::vector,y::vector) 
local i,j,k,n; global T, alpha, beta, g,z,X; 

n:=vcctdim(d):alpha:=array(l..n-l):beta:=array(l..n): g:=array(l..n): 
z:=array(l..n):X:=array(l..n): 

^components of the vectors alpha, beta, and gamma # 
beta[l]:=d[n]: 

if beta[l]=0 then beta[l]:=x;d[n]:=x;fi: 

g[2]:=a[n-l]/beta[l]:alpha[l]:=b[n]: beta[2]:=simpUfy(d[n-l]-alpha[l]*g[2]): 
alpha[2]:=simpUfy(b[n-l]-g[2]*bb[n]): 
if beta[2]=0 then beta[2]:=x;fi: 
for i from 3 to n-1 do 

g[i]:=simplify((a[n-i+l]-aa[n-i+l]*alpha[i-2]/beta[i-2])/beta[i-l]): 

alpha[i] :=simplify (b[n-i+ l]-g[i] *bb[n-i+2] ) : 

beta[i]:=simplify(d[n-i+l]-aa[n-i+l]*bb[n-i+3]/beta[i-2]-alpha[i-l]*g[i]): 
if beta[i]=0 then beta[i]:=x;fi: 

end do: 

g[n]:=simplify((a[l]-aa[l]*alpha[n-2]/beta[n-2])/beta[n-l]): 
beta[n]:=simplify(d[l]-aa[l]*bb[3]/beta[n-2]-alpha[n-l]*g[n]): 

# To compute the vector Z # 
z[l]:=y[n]:z[2]:=y[n-l]-g[2]*z[l]:i:='i': 
for i from 2 to n do 

z[i]:=simplify(y[n-i+l]-aa[n-i+l]*z[i-2]/beta[i-2]-g[i]*z[i-l]): 

end do: 

To compute the Solution of the system X. 
X[n]:=z[n]/beta[n]:i:='i': 

X[n-l]:=simplify((z[n-l]-alpha[n-l]*X[n])/beta[n-l]): 
for i from n-1 by -1 to 1 do 

X[i]':=simplify((z[i]-alpha[i]*X[i+l]-bb[n-i+l]*X[i+2])/beta[i]): 

end do: 

# To compute the determinant T # 
T:=subs(x=0, simplify (product (beta[r],r=l..n))); 
eval(X): 

end: 

> # Call no. 1 for the procedure back_penta. # 

> x:='x': 

> aa:=vcctor([3,-l,7,-2]); 

aa := [3, -1, 7, -2] 
>a:=vcctor([2,5, 2,3,-5]); 

a := [2, 5, 2, 3, -5] 

> d:=vcctor([l,3,3,5,6,14]); 

d := [1, 3, 3, 5, 6, 14] 

> b:=vector([0,2,l,2,2,l]); 

b := [0, 2, 1, 2, 2, 1] 

> bb:=vector([0,0,-5,-7,3,-10]); 

bb := [0, 0, -5, -7, 3, -10] 

> y:=vector([6,9,8,l,6,5]); 

y := [6, 9, 8, 1, 6, 5] 

> X:=back_penta(bb,b,d,a,aa,y); 

^ X:= [1, 1, 1, 1, 1, 1] 

>T; 

-8597 
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> # End of call no. 1. # 



># Call no. 2 for the procedure back_tri. # 

> x:='x': 

> aa:=vector( [3, -1,7,-2]); 

aa := [3, -1, 7, -2] 
>a:=vector([2,5, 2,3,-5]); 

a := [2, 5, 2, 3, -5] 

> d:=vcctor([l,3,3, 5,6,0]); 

d := [1, 3, 3, 5, 6, 0] 

> b:=vcctor([0, 2, 1,2,2,1]); 

b [0, 2, 1, 2, 2, 1] 

> bb:=vector([0,0,-5,-7,3,-10]); 

bb := [0, 0, -5, -7, 3, -10] 

> y:=vector([6,9,8,l,6,-9]); 

y := [6, 9, 8, 1, 6, -9] 

> X:=back_penta(bb,b,d,a,aa,y); 

v r 1777 2z+1777 489a;-1777 116(to-1777 574a-1777 182^ + 1777 ] 

A - — L 741x-1777' 741rr-1777' 741a;-1777' 741a;-1777 ' 741^-1777' 741x-1777J 

> T; 

1777 

> x:=0:X:=map(eval,op(X)); 

X:=[l, 1,1, 1,1,1] 

> # End of call no. 2. # 
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